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ABSTRACT 

In order to know the formation epoch of the oldest elliptical galaxies as a 
function of mass and observed redshift, a statistical analysis for 333 extremely 
red objects (EROs) classified as old galaxies (OGs) at 0.8 < z < 2.3 is carried 
out. Once we get My and (B—V) at rest for each galaxy, we calculate the average 
variation of this intrinsic color with redshift and derive the average age through 
a synthesis model (the code for the calculation of the age has been made publicly 
available). The average gradient of the (B — V) color at rest of EROs/OGs is 
0.07-0.10 Gyr -1 for a fixed luminosity. The stars in these extremely red elliptical 
galaxies were formed when the Universe was ~ 2 Gyr old on average. We have not 
found a significant enough dependence on the observed redshift and stellar mass: 
f ormation = _q A q ± ( ^Mfe = _o.81 ± 0.98) Gyr. This fits a scenario 

\ ("observed /' \dlog w M, J J 

in which the stellar formation of the objects that we denominate as EROs-OGs 
is more intense at higher redshifts, at which the stellar populations of the most 
massive galaxies form earlier than or at the same time as less massive galaxies. 

Subject headings: galaxies: formation — galaxies: high redshift — galaxies: 
statistics — Infrared: galaxies 



1. Introduction 

It is usually sustained that galaxies at high redshift are intrinsically bluer than at low 
redshift (e.g., Dickinson et al. 2003; Rudnick et al. 2006; Labbe et al. 2007), as would 
expected if their populations were younger and with lower mass/ luminosity ratios. However, 



-2- 



analysis of this intrinsic color evolution is not free from caveats, due mainly to the difficulty 
in disentangling selection effects. For instance, Dickinson et al. (2003, Fig. 2) showed a 
clear lack of red objects between 2 < z < 3 with ( m 1700 ^ ~~ m B)AB,rest > 3, while there were 
many of these red objects at z < 2; so they consequently claimed that there is an evolution 
in color. However some red galaxies were probably missed in that analysis. As a matter of 
fact, if we take the subsample of old elliptical galaxies without extinction (/ id > 0.95) with 
2 < z < 3 from Miyazaki et al. (2003) (N = 12 galaxies) we find that all of them have 
2.5 < ( m 170a | — mB)rest < 4.5 with an average of 3.6. Miyazaki et al.'s galaxies are much 
redder than Dickinson et al.'s, which shows that the Dickinson et al. (2003) sample does not 
contain the reddest galaxies. 

On the question of the ages of galaxies and when were they formed, there also some 
uncertainties. There are many observations at different redshifts of galaxies nearly as old 
as the Universe, although the density of those galaxies at high redshift is significantly lower 
(Renzini 2006; Abraham et al. 2007). Early-type massive galaxies with early formation were 
found at z ~ 1 (di Serego Alighieri et al. 2006; 2007), z = 1 — 2 (Spinrad et al. 1997; Daddi 
et al. 2005; Longhetti et al. 2005; Toft et al. 2005; Trujillo et al. 2006), z = 2-3 (Daddi 
et al. 2005; Labbe et al. 2005; Cassata et al. 2008, Kriek et al. 2008, 2009), z = 3 - 4 
(Toft et al. 2005; Chen & Marzke 2004), z = 4 - 5 (Chen & Marzke 2004; Rodighiero 
et al. 2007) and z > 5 (Wiklind et al. 2008). Most of this information is obtained from 
photometry, but there are also some spectra massive elliptical galaxies at z — 1.4 — 2.2 
(Cimatti et al. 2008, Kriek et al. 2009) revealing them to be old galaxies. Semi-analytical 
ACDM models (De Lucia et al. 2006) claim that the formation of low mass galaxies is first to 
give way to mergers of very massive galaxies, which is apparently at odds with observation 
of these massive galaxies at high redshift, even if they were formed through dry mergers in a 
downsizing scenario. Ferreras et al. (2009) also found that very massive galaxies do not have 
significant evolution at z < 1.2. However, Schiavon et al. (2006), who have taken spectra of 
red galaxies (U — B > 0.25) at z ~ 0.9, derived their age to be on average 1.2 Gyr, which 
means that some galaxies were formed at lower redshifts. Arnouts et al. (2007) also show 
evidence for a major build up of the red sequence between z = 2 and z — 1. 

In this paper, we pay further attention to the determination of intrinsic color and age 
variation for different redshifts in a statistical way for two different samples of very red 
elliptical galaxies in order to constraint their formation epoch. 
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2. Data 

We select galaxies classified as Extremely Red Objects (EROs) within the redshift range 
0.8 < z < 2.3 and, within this group, those classified as Old Galaxies (OGs) with negligible 
intrinsic extinction; that is, passively evolving populations of elliptical galaxies. The different 
methods of selecting OGs are quite consistent with each other (Fang et al. 2009). Galaxies 
will be selected with available fluxes in the three near infrared filters (JHK), plus at least 
two filters in the optical for which the flux signal/noise is greater than 3. We use two sources 
of publicly available data: 

1. ECDFS catalog (Taylor et al. 2009a), which gives photometry in ten filters: U, U38, B, 
V, R, I, z', J, H, K s from ISAAC (VLT)+Hubble data. For the selection of EROs-OGs, 
we adopt the color criterion (i 775 — K)ab > 2.42 (EROs), (J — K)ab < 0. 20(2775 — 
K)ab + 0.39 (Fang et al. 2009). We do not have the magnitude at 2775 but we get 
this with the corresponding color correction using the adjacent filters. A total of 276 
galaxies. 

2. Miyazaki et al. (2003) give photometry in eight filters: B, V, R, i', z', J, H, K s 
from Subaru/XMM-Newton+UH2.2m. Their sample of EROs-OGs was selected using 
(R — K)ab > 3.35 (EROs), and within these sources, by means of spectrum fitting, 
as OGs without extinction with a fraction of old population j Q u < 0.95. A total of 57 
galaxies. 

Estimating roughly the average intrinsic color/age of EROs/OGs is our aim here. The 
term "ERO" reflects the observed characteristic color of a galaxy, not its intrinsic properties. 
For this reason, and because of using magnitude-limited samples, they have different ranges of 
stellar masses, M/L ratios and intrinsic colors at different redshifts. It is not a homogeneous 
sample of galaxies with the same intrinsic characteristics at all redshifts; there are biases. 
Nonetheless, throughout the paper we shall separate the dependence on redshift from the 
dependence on luminosity/mass. 

3. Color (B - V) at rest 

We take AB apparent magnitudes, corrected for Galactic extinction (although negligi- 
ble), for different wavelengths, m^B^), (i = Aj), for 5 < Nf < 10, with the corre- 
sponding error bars (in our case in optical and near infrared). As has already been said, we 
consider only the points with a flux signal/noise above 3. We will use data with available 
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redshifts, most of them photometric. The average systematic error of the photometric red- 
shifts is 0.025 for ECDFS (Taylor et al. 2009a, §7.3), which is small, so we do not 

take it into account here; similarly for the Miyazaki sample. There is a statistical error for 
each photometric redshift, but we expect that these uncertainties will nearly cancel in the 
statistical analysis. 

With this, we calculate the rest luminosities in two filters. This is done through spectral 
energy distribution (SED) fitting using templates of galaxies with the software Interrest v2.0 
(Taylor et al. 2009a). The calculations for the sample ECDFS have already been carried 
out by Taylor et al. (2009a). The calculation for Miyazaki et al. (2003) was carried out 
by us with the Interrest v2.0 software, by changing to the Subaru filters. In this paper, 
we do the calculations only with a pair of filters at rest (Johnson B and V), applying the 
correction to convert AB into the Vega calibration: (B — V) = (B — V)ab + 0.119 (Frei 
& Gunn 1994). The (U — B) rest color is not used in this paper because it is more sensitive 
to redshift uncertainties and uncertainties in the emission-line corrections (Rudnick et al. 
2006). Moreover, the (U — B) color depends more strongly than (B — V) on metallicity for 
a given age; it is also more sensitive to a-enhancement (Cassisi et al. 2004). 

There are uncertainties in the rest color, due to the error bars of the apparent mag- 
nitudes, the deviation of the assumed shape from the true SED (spectral features which 
move through the filter bands), errors in the photometric redshifts, etc. In any case, we do 
not expect important systematic errors, and the statistical errors can be reduced when we 
calculate the average for bins with a large number of galaxies. 

In Fig. [TJ we give the results of the average colors as a function of the age of the Universe 
when the galaxy is observed, 



with Ho = 73 km/s/Mpc, Q m = 0.24, = 0.76. Both samples give approximately the same 
results; thus the selection of OGs among EROs is shown to be quite consistent with both 
independent methods. There is a significant average gradient in color of 0.124 ± 0.014 and 
0.118 ± 0.021 Gyr -1 respectively for samples ECDFS and Miyazaki et al. 

Note that the galaxy gap in the lower right corner of Fig. [T] is at least partially an 
artifact of the sample intrinsic color bias as a function of the redshift. This is expected 
because at the lowest redshift, the (i — K) and (R — K) limits used will pick out only the 
very oldest and reddest galaxies, whereas, at higher redshifts, younger and bluer galaxies 
will be included. Note even that at the very lowest redshifts in Fig. [1] the much smaller 
Miyazaki sample dominates; this is because the ECDFS criterion is a bit more strigent, so 
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it eliminates even the oldest passive galaxies at z ~ 0.8. One must therefore interpret Fig. 
[l]as the intrinsic color as a function of redshift of the objects selected as EROs/OGs, not of 
a general characterization for elliptical galaxies. 

If we do a bi-linear fit of the colors as a function of two independent variables t b s . and 
M v>rest , we get: 

(B - V) rest = ai + a 2 [t o6s .(Gyr) - 5] + a z {M v>rest + 22), (2) 

with ai = 0.756±0.014, a 2 = 0.084±0.015, a 3 = 0.065±0.013 for ECDFS; a x = 0.814±0.023, 
a 2 = 0.087 ± 0.027, a 3 = 0.048 ± 0.028 for Miyazaki. This bi-linear fitting allows us to 
separate the evolution from the biases in absolute magnitudes; so the second term (0.084 
or 0.087) gives us the average evolution in color for a fixed luminosity. Galaxies are redder 
for lower observed redshift and for lower luminosity. The first fact indicates older galaxies 
at lower redshift, and the second fact is probably related to a higher luminosity for younger 
populations. Other authors, for instance Labbe et al. (2007) have found however that the 
most luminous galaxies have redder colors, with a slope of d ^^^ = —0.09 ± 0.01 (Labbe 
et al. 2007). Our guess is that we do not find the same luminosity dependence because we 
have preselected massive galaxies with the constraint that they be EROs/OGs, and within 
a different range of redshifts. In Taylor et al. (2009b, Fig. 4), we see that d ^~ r - > is negative 
for z < 1.25; however, the trend changes for z > 1.25 and this variation of average color with 
absolute magnitude is null or even slightly positive. 



4. Age estimation for early-type galaxies 

NOTE: A FORTRAN code to carry out the calculations explained in this section is 
available at \http:/ / www, iac. es/ galeria/martinlc/ codes. htm\ 

As has been said, we have selected old massive elliptical galaxies with negligible internal 
extinction, i.e. without gas and dust. We can therefore connect the color and luminosity 
in V-rest of the galaxy with the average age of the stellar population and its metallicity. 
There may be some wrong identification of OGs among our galaxies selected with the color 
method for the ECDFS sample (Miyazaki et al. 2003 and Fang et al. 2009 estimate it to 
be ~ 25%), and consequently there may be some case of dusty galaxies among our ECDFS 
sample, but the statistical comparison with the SED fitting method of the Miyazaki et al. 
sample in Fig. [T] shows that there are neither significant differences nor systematic effects 
with redshift. Only perhaps in the range 1.2 < z < 1.6 might there be some small difference, 
where the starburst contamination might be higher (Fang et al. 2009). 

In order to estimate the average age corresponding to our galaxies, we use a synthesis 
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model: Vazdekis et al. (1996; hereafter V96); see also Vazdekis (1999). There is an age- 
metallicity degeneracy, but this can be broken approximately with the use of the mass- 
metallicity correlation. Another way to break the degeneracy would be by using two colors 
(Li & Han 2007, and references therein), but we have only one reliable color at rest and we 
do not have rest near-infrared colors as necessary in Li & Han (2007). 

We must also bear in mind that synthesis models return the mean value of a distribution, 
and a perfect fit to observational data to infer the age is only correct on average, since 
the individual cases may present some dispersion with respect to the average (Cervino & 
Luridiana 2006). For this reason, and because of the large errors in the color of each galaxy, 
we do not calculate the age of each galaxy separately but the average age of each bin of 
galaxies (11 in our case) with the same redshift. 

For each bin of galaxies, the steps are as follows: 

1. We assume zero metallicity and derive the age (t\) of the galaxy which is given by the 
V96 model for the given (B — V) TCSt of the galaxy. 

2. Given the age t\ and the zero metallicity, we derive with the V96 model the stellar 
mass-to-light ratio in the V filter [(M*/Ly)i]. 

3. Since we know the luminosity at rest in the V filter (Ly), we can derive the stellar 
mass of the galaxy: 

M„,i = (M./Lv^Lv. (3) 

4. Given the stellar mass of the galaxy, we estimate the metallicity [Fe/H]i. We use the 
correlations of metallicity and a-enhancement given by Thomas et al. (2005). The 
average relation is: 

[Fe/H] = 0.066 log in ( ) ±0.11, (4) 

with and uncertainty of 8 ([Fe/H]) ~ 0.11 including the scatter of the correlations 
and the variations with the environment (low or high densities) of the galaxies. The 
correlation of mass (or velocity dispersion) with metallicity is also observed in Yamada 
et al. (2007). There is also a dependence on age, but it is only important for low mass 
galaxies of velocity dispersion less than 100 km/s (Yamada et al. 2007), which is not 
the case of our galaxies. It is assumed that the relationship of mass and metallicity 
does not evolve with redshift in a passive evolution (di Serego Alighieri et al. 2006). 
A metallicity evolution in galaxies of the same mass is not observed (Cimatti et al. 
2008) so, provided that mass does not correlate with the age, the mass-metallicity 
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relationship should not change too much at high redshift. As is observed in Fig. [2], the 
accuracy in the metallicity determination mainly affects the reddest (oldest) galaxies. 
For the youngest galaxies (< 3 Gyr) the errors in metallicity are not very important, 
so a possible evolution in the relationship of Eq. (J3J) at high redshift would not affect 
the results. 

5. We derive again the age t 2 with the color (B — V) rest and metallicity [Fe/H]i. 

6. We repeat steps 2-5, which gives an age £3. Since £3 ~ t 2 , we do not need to do further 
iterations and we have got the convergence of the metallicity, mass and age with only 
three iterations. If £3 were significantly different from t 2) we would continue to iterate 
to get £4, £5,... until the convergence is obtained. 



There is some dependence with the IMF (initial mass function) slope. In Fig. [2l we 
plot age vs. color for metallicities [Fe/H] = and [Fe/H] = 0.2 and different slopes in the 
bimodal IMF, as defined in V96, or Kroupa IMF using V96 model. Slope 1.3 in the bimodal 
IMF is the standard value and is nearly coincident with the Kroupa case. Within variations 
of ±1 of the slope the variations of the age are fitted from Fig. [2] for null metallicity by: 

i^S— = -0.071 + 0.346(5 - V) - 0.257(5 - V) 2 , 

^j- = 0.123 - 0.299(5 - V) + 0.243(5 - V) 2 . (5) 

These variations produce some error in the age determination. However, this error is smaller 
than that produced by the uncertainty in the color or the metallicity. In any case, they are 
taken into account. 

The age of the early-type galaxies is plotted in Fig. [31 The vertical bars include the 
errors due to uncertainties in the colors, the uncertainty of 0.11 in the mass-metallicity 
relationship of Eq. (T4J) and the variations due to the IMF slope change within a range of ±1 
given by Eq. ([5]). There is no zero-point calibration problem except perhaps for the last bin 
because this affects mainly galaxies older than 5 Gyr (Vazdekis et al. 2001). We must also 
bear in mind that we are neglecting the systematic errors in the photometric redshifts; were 
they non-negligible, we would have extra systematic errors in the calculated ages. 

Figure [3] represents the average age of the given sample among the ERO/OGs, with 
all the selection effects associated with each redshift. If we separate the dependence on 
mass from its evolution, subdividing each redshift bin into sub-bins with different luminosity 
(AMy = 1), and we do over them a bilinear fit weighted with the square inverse of relative 
errors, both for average color and average age, we get 
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Fig. 1. — (B — V) Test (Vega calibration) colors observed at different ages t b s .- 
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Fig. 2. — Variation of age with (B — V) color in the V96 stellar population synthesis model. 
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(B - V) rest = h + b 2 (t ohs . - 5) + 6 3 k>g 10 (i\4) , (6) 

igai. = ci + c 2 (t ohs . - 5) + c 3 log 10 (Af«) , (7) 

with t ga i and tf orm . in units of Gyr, M* in units of 10 11 M Q ; b% = 0.779 ± 0.025, b 2 = 
0.106 ± 0.014, 63 = -0.042 ± 0.049, a = 2.22 ± 0.54, c 2 = 1.15 ± 0.32, c 3 = -0.4 ± 1.0 for 
ECDFS; and h = 0.845 ± 0.020, b 2 = 0.085 ± 0.022, 63 = 0.0027 ± 0.0020, a = 5.06 ± 0.96, 
c 2 = 2.71±0.78, C3 = 4.0±2.1 for Miyazaki. Average stellar masses in the bins range between 
5 x 10 9 Mq and 3 x 10 11 M . Kriek et al. (2008) found in the same range of redshifts a 
d( - u ~ B ^ Icst = +0.024 Gyr" 1 for the red sequence at a given mass of 2 x 10 11 M Q , equivalent 

^obs. 

(with « 0.8; V96 at [Fe/H]=0, {U - B) = 0.25-0.37) to b 2 = « +0.03 Gyr" 1 , 

a smaller color evolution but for larger masses than our sample. 



5. Discussion on their formation epoch 

The average epoch of star formation (the epoch of formation of the first stars, might 
be lower) is tf orm _ = t Q b s . — t ga \_, shown in Fig. HI Separating the evolution from the mass 
dependence, 

korm. = d l + d 2 (t ohs , -5) + d 3 log 10 (M*), (8) 

withdi = 1.94+0.51, d 2 = -0.46±0.32, d 3 = -0.81±0.98 on average for the ECDFS+Miyazaki 
samples. The observed EROs/OGs are all formed within a narrow range of epochs, when 
the Universe was less than 4 Gyr old (z > 1.7). The present fit is compatible with all 
EROs/OGs formed at the same time at, on average, tf or m. = 2.0 ± 0.3 Gyr (z ~ 3 — 4). 
Again, we remind the reader that the criterion to select EROs/OGs is more restrictive at 
lower redshift, picking out only very old galaxies, while at higher redshift the range of al- 
lowed ages is wider. We might therefore expect that the oldest EROs at high z (low t b s .) 
will have an even lower formation age. The age tf or m. = 2.0 Gyr is a conservative lower limit 
representing the average sample; there must be some EROs/OGs formed beforehand. Given 
that ^j*™- = -0.46 ± 0.32 « +1 for a given stellar mass, this means that the galaxies of 
the present sample are not formed continuously with the same rate, but more intensely at 
higher redshifts. 

The stellar populations of most massive galaxies are not formed much later than the 
less massive ones (d 3 < 0). This agrees the results mentioned in the introduction that very 
massive evolved galaxies detected at redshifts 1.5-6 were formed in the very early Universe 
(Daddi et al. 2005: Chen & Marzke 2004; Rodighiero et al. 2007; Wiklind et al. 2008). 
This might appear in contradiction with the result of §3] that galaxies are redder for lower 
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luminosities, but it does not. As said, the mass-luminosity ratio does not remain constant 
giving higher luminosity for younger objects so it is not contradictory that older/redder 
objects correspond to lower luminosities and higher masses. It is in fact observed if we 
compare Figs. 4 and 5 of Taylor et al. (2009b) for z > 1.25: clearly, a strong dependence on 
stellar mass does not mean a strong dependence on luminosity. 
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